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-■»«. ■ Abstract 

Ph, Bulk properties of strange quark matter (SQM) are investigated within the SU(3) 

Q_i| Nambu-Jona-Lasinio model. In the chiral limit the model behaves very similarly 

to the MIT bag model which is often used to describe SQM. However, when we 
introduce realistic current quark masses, the strange quark becomes strongly disfa- 
vored, because of its large dynamical mass. We conclude that SQM is not absolutely 
^ ■ stable. 
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The properties of strange quark matter (SQM) and in particular the conjecture that 
SQM could be the absolute ground state of strongly interacting matter |Lj, || have at- 
tracted much attention in nuclear physics and astrophysics. Whereas the empirical sta- 
bility of ordinary nuclei excludes the existence of absolutely stable non-strange quark 
matter (NSQM) this argument does not hold for SQM because the decay of nuclei into 
SQM would involve higher-order weak interaction processes, associated with a very long 
lifetime. In 1984 Farhi and Jaffe |§ investigated the question of absolutely stable SQM 
within the MIT bag model [|J. Treating the bag constant and the strange quark mass as 
free parameters they found a reasonable window for which SQM is stable and NSQM is 
unstable compared with an 56 Fe-nucleus. Almost 15 years later bag model calculations 
still play a key role in this field (for reviews see [f|-|§). 

In ref. we have investigated non-strange quark matter within the Nambu-Jona- 
Lasinio (NJL) model |HJ in mean field approximation. For vanishing (current) quark 



masses and with the appropriate choice of parameters we found that quark matter in the 
NJL mean field behaves very similarly to quark matter in an MIT bag. For the description 
of SQM we have to introduce finite quark masses, at least in the strange quark sector. 
In this case the correspondence between the two models becomes less accurate, with the 
NJL mean field behaving in a much more complex way. This is related to the fact that 
the bag constant in the NJL model is not a phenomenological input parameter, like in a 
bag model, but a dynamical property, which has its origin in the spontaneous breaking 
of chiral symmetry. This motivated us to extend the model of ref. |J to investigate bulk 
properties of strange quark matter within the NJL model. 

The three-flavor NJL model has been discussed by many authors, e.g. [ITTJ- |PT[ . In 



this article we adopt the Lagrangian of ref. [21 



£ = q(i$-m )q + G^[{q\ k qf + {qi^\ k qf} 

- K[det f (q{l + l5 )q)+det f (q(l- l5 )q)} . (1) 

Here q denotes a quark field with three flavors, u, d and s, and three colors. rfiQ = 
diag(mQ,rriQ,rn S Q) is a 3 x 3 matrix in flavor space. We restrict ourselves to the isospin- 
symmetric case, ttIq = m^. The Lagrangian contains a U(3)l x £/(3)^-symmetric four- 
point interaction term and a six-point interaction, which is a determinant in flavor space 
and which breaks the £7A(l)-symmetry. G and K are constants with dimension energy -2 
and energy -5 , respectively. 

The model is not renormalizeable and we have to specify a regularization scheme for 
divergent integrals. For simplicity we use a sharp cut-off A in 3-momentum space. Besides 
A we have to fix four parameters, namely the coupling constants G and K and the current 
quark masses mj and m s Q . In most of our calculations we will adopt the parameters of 
ref. 0: A = 602.3 MeV, GA 2 = 1.835, KA 5 = 12.36, m% = 5.5 MeV and m s = 140.7 MeV. 
These parameters have been determined by fitting m n , mx and m v > to their empirical 
values, while the mass of the f^-meson is underestimated by about 6%. 

Due to the interactions the quarks acquire dynamical masses rrti which are in gen- 
eral different from their current masses m l and depend on temperature and density. In 
the following we consider quark matter at T = and zero or non-zero baryon number 
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density ps = \ub = \{n u + nd + n s ), with n« = (qjqi)- Then, in mean field (Hartree) 
approximation, the dynamical quark masses are solutions of the gap equation 

rrii = m - 4G (q^i) + 2K (q j q j ){q k q k ) , (2) 

with qj, q k ) being any permutation of (u, d, s). The quark condensates {q%qi) are given 

by 3 A 

(QiQi) = 2 P 2(i P i „ 1 ■ ( 3 ) 

They depend on the Fermi momenta Pf% = (tt 2 ^) 1 ^ 3 an d again on the dynamical masses 
mj. Thus, for given quark number densities rii, eqs. (0) and ([|) have to be solved self- 
consistently. After that we can calculate the energy density e and the total pressure p of 
the system: 



3 fPFi j I 

e= — / ' P 2 dp Jmf +p 2 - (B-B ) , p=-e+ ^ ^V^f+PFf- (4) 



3 /"Pf 

Here i? denotes the bag pressure 
3 ' A 

{ 



B = E ( -5 / ^ (\/ m i+P 2 ~ ^o+P 2 ) ~ 2G(q t , qi } 2 ) + 4K(M M )(Jrf)(s S ) , 

(5) 

and So = -B|n u =n d =n s =o- The latter was introduced in eq. (0J) in order to ensure e = p = 
in vacuum. 

The existence of a bag pressure term in eq. (0) suggests that there is a certain con- 
nection between the present model and standard bag model descriptions of quark matter. 
Let us have a closer look at this point. It is useful to begin with the chiral limit, i.e. 
m = for all flavors. With A, G and K as specified above, chiral symmetry is spon- 
taneously broken in vacuum with equal dynamical masses m u = ma = m s = 310.6 MeV 
and a bag pressure Bq = 57.3MeV/fm 3 . (There are other solutions of the gap equation 
but they are energetically less favored.) Now we consider S*?7(3)-symmetric quark matter, 
n u = n d = n s = Pb- At low densities the dynamical quark masses decrease with density 
but they remain finite. However, when ps exceeds a critical value of ~ 0.9po chiral sym- 
metry gets restored and the quark masses and condensates vanish. Consequently the bag 
pressure is also zero in this regime and we get from eq. (^) 

9 3 
-p F A + B , p = —^Pf 4 - Bq . (6) 



This means the system behaves exactly like a gas of massless non-interacting quarks with 
three flavors and three colors inside a large MIT bag with a bag constant Bq. It should 
be kept in mind, however, that eq. @ only holds in the regime of vanishing dynamical 
quark masses. For the present parameters it turns out that the minimal energy per baryon 
number E/A = e/ps is indeed found inside this region. Thus eq. (H) is appropriate to 
describe the system in the vicinity of the point of stability. 

For comparison we consider non-strange isospin-symmetric quark matter, n s = 0, 
n u — n d — \pb-i still in the chiral limit. In this case strange and non-strange quarks 
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behave differently with density, and chiral symmetry is only partially restored at high 
densities: When p B exceeds ~ 1.0p the non-strange quarks become again massless but 
the strange quark remains massive with a mass m* which is 154.2 MeV in our example. 
Hence the bag pressure remains also finite, B = B* = 3.7MeV/fm 3 , and the system now 
behaves like a bag model for a gas of massless quarks with two flavors and three colors 
but with an effective bag constant B e ff = B — B*, which is smaller than the bag constant 
Bq in the SU(3)-symmetric case. 

Thus, in contrast to a naive bag model, the effective bag constant which describes the 
behavior of NJL quark matter depends on the flavor composition of the matter. This 
is also shown in fig. 1 where various properties of isospin-symmetric quark matter are 
plotted as functions of the fraction of strange quarks, r s = n s /n B . With the present 
parameters stable quark matter exists in the region < r s < 0.48. As discussed in ref. || 
these solutions correspond to the high-density phase of a first-order chiral phase transition 
in equilibrium with the non-trivial vacuum. Throughout this paper we will denote the 
dynamical quark masses of these stable matter solutions by m* and the corresponding 
bag pressure by B*. In fig. 1(a) the solid line indicates the energy per baryon number, 
E/A of stable quark matter, m* and the effective bag constant -B e // = Bq — B* are shown 
in fig. 1(b) and (c), respectively. For 0.12 < r s < 0.48 the system behaves like a bag 
model with massless quarks and a bag constant -B e // = B . For r s < 0.12 the effective 
bag constant is somewhat smaller and the strange quarks become massive. In panel (a) 
we also show the bag model result for E/A with a bag constant B BM = B and massless 
quarks (dashed line). For 0.12 < r s < 0.48 it agrees of course exactly with the NJL result, 
but also for r s < 0.12 the deviations are relatively small. In fig. 1 all NJL curves stop at 
r s ~ 0.48. For r s > 0.48 the NJL quark matter is unstable against evaporation of free 
massive s-quarks. In fact, this kind of phase separation is very typical for asymmetric 
matter and has been discussed for nuclear matter [^2], [23j as well as for models of the 
deconfinement phase transition |24|, p5[| . However, the evaporation of free quarks is of 
course unphysical and reflects the missing confinement of the NJL model. 

With finite current quark masses m l Q things change considerably. Fig. 2 shows the 
same quantities as fig. 1, but now with tjIq = tjIq = 5.5 MeV and m s = 140.7 MeV [21]. 



We now find stable quark matter for < r s < 0.79. In vacuum the dynamical masses of 
the up- and down-quarks are increased by ~ 60 MeV to 367.6 MeV and that of the strange 
quark by ~ 240 MeV to 549.5 MeV. Perhaps more important for us is the well-known 
fact that the quark condensates do not completely vanish even at higher densities. As a 
consequence, the dynamical quark masses in stable quark matter, m*, are still well above 
the current masses m l (see fig. 2(b)). For non-strange isospin-symmetric quark matter 
(r s = 0) we find m* = 52.6 MeV and m* = 464.4 MeV. Since on the other hand the Fermi 
energy in this system is only 361.1 MeV the conversion of a non-strange quark into a 
strange quark is obviously energetically not favored. Here the notion of dynamical quark 
masses is crucial: In a naive bag model one would assign the strange quark its current 
mass which in our example is more than 200 MeV lower than the Fermi energy. 

The difference becomes obvious in fig. 2(a) where the energy per baryon number of 
stable quark matter is plotted as a function of r s . The solid line indicates the NJL 
result while the dotted line corresponds to a bag model calculation with B B m = 105.2 
MeV/fm 3 (the effective bag constant at r s = 0) and quarks with m u = = 5.5 MeV 
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and m s = 140.7 MeV, our current quark masses. Whereas the latter has a minimum at 
a finite fraction of strange quarks, r s = 0.27, the former is a strictly rising function of 
r s : additional strangeness is always disfavored in the NJL calculation. In fig. 2(a) we 
also show the result of a bag model calculation with m u = rrid = 52.6 MeV and m s = 
464.4 MeV, our dynamical quark masses at r s = (dashed-dotted line). This curve rises 
even steeper with r s than the NJL result. The reason is that in the NJL calculation the 
dynamical strange quark mass m* drops with r s (see fig. 2(b)). Therefore the penalty for 
adding strangeness becomes smaller than it would be with a constant mass m*(r s = 0), 
even though this effect is partially compensated by the fact that the effective bag constant 
B e ff = Bo — B* rises with r s for not too large r s . 

B e ff is shown in fig. 2(c). It varies much stronger with r s than in the chiral limit 
and is much larger in magnitude. Already at its minimum, at r s = 0, it is about twice 
as large. To major extent this is the reason why the energy per baryon number at this 
point is almost 200 MeV larger than in the chiral limit. Both, the magnitude of the 
effective bag constant and its strong variation with r s , have their origin in the fact that 
the bag pressure is very sensitive to the dynamical quark masses. As a consequence the 
vacuum bag pressure B is now 291.7 MeV/fm 3 , five times the chiral limit value, and, 
since the strange quark mass in stable matter varies, the corresponding bag pressure B* 
varies strongly with r s . 

It should be also noted that, even for fixed r s , quark masses and the bag pressure 
are density dependent and the values for m* and B e ff shown in fig. 2 only correspond to 
the densities of stable quark matter. As soon as one moves away from the minimum of 
E/A, rrii and Bo — B change. Thus, whereas in the chiral limit NJL quark matter could 
be well described in terms of an MIT-like bag model with vanishing quark masses and a 
single bag constant (at least within a larger range of densities and strangeness fractions), 
it behaves in a much more complex way when we include finite current quark masses. 

So far, focusing on the bag model aspect, we made some simplifying assumptions 
which we have to abandon for a more realistic study of SQM. In particular we have to 
take into account weak decays. This implies that we have to include electrons and (in 
principle) neutrinos. To large extent we will adopt the model of Farhi and Jaffe ||, 
with the MIT bag replaced by the NJL mean field. With the usual assumption that the 
neutrinos can freely leave the system, the matter is characterized by four (three quark- 
and one electron-) chemical potentials. They are related to the corresponding densities 
in the standard way: 



rii = — (n? — m\ ) 3//2 9(/j^ — ml) for i = u,d,s and n e = . (7) 



Here we neglected the electron mass. For the quarks one has to keep in mind that the 
masses entering the r.h.s. have to be calculated from the gap equation and are therefore 
density dependent themselves. Treating the electrons as a free gas, energy density and 
pressure of SQM are given by 



7T 



1 
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with e and p as defined in eq. 
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In chemical equilibrium maintained by weak interactions only two of the four chemical 
potentials are independent: 

= Us = Hu + Pe • (9) 
Furthermore we demand charge neutrality, 

2 1.. 

^n u - -(n d + n s ) - n e = . (10) 

Thus the system can be characterized by one independent quantity, e.g. the baryon 
number density p#. 

Our results are shown in fig. 3. Panel (a) shows the energy per baryon number E/A = 
£sqm/pb as a function of p#. The solid line corresponds to the model described above, 
the dotted line to non-strange quark matter, where p s was set equal to zero by hand. 
Obviously the strangeness degree of freedom is only important at densities above ~ 4 p . 
This can also be seen in panel (c) where the fractions = rij/rif, of the various particles 
are plotted. Since electrons (dotted line) play practically no role (and are hardly visible in 
the plot) the fraction of u-quarks (dashed-dotted) is fixed by charge neutrality (eq. fllPP) 
to r u ~ 1/3. For p < 3.85 p the fraction of strange quarks, r s , (solid) is zero and thus r d 
must be the remaining 2/3 (dashed). For p > 3.85 po r s becomes non-zero and r d drops 
accordingly. 

The fact that there is no strangeness at lower densities is again due to the relatively 
large mass of the strange quark. The dynamical quark masses are plotted in fig. 3(b) as 
functions of ps- For comparison we also show the chemical potential p s = p d (dotted 
line). As long as p s < m s the density of strange quarks is zero. For p B < 2p , m s drops 
with density. This can be mainly attributed to the decrease of {uu) (dd) due to the rising 
density of non-strange quarks. At higher densities (uu)(dd) is practically zero and m s 
stays almost constant until n s becomes non-zero and causes \{ss) \ to drop. This behavior 
is rather different from standard parameterizations which have been used in the literature 
to study SQM and which depend on the total baryon number density ps only |26], 27 \. 

At those densities where strange quarks exist, they lead to a reduction of the energy, 
as can be seen by comparing the solid curve in fig. 3(a) with the dotted one. However, the 
minimum of E/A does not lie in this regime, but at a much lower density, ps = 2.25 po- 
Here we find E/A = 1102 MeV. Compared with the energy per baryon in an iron nucleus, 
E/A ~ 930 MeV, this is still very large. In this sense our results are consistent with the 
empirical fact that stable NSQM does not exist. However, since the energy of strange 
quark matter is even higher, our calculation predicts that also SQM is not the absolute 
ground state of strongly interacting matter. 

This result is very robust with respect to changes of the model parameters. Obviously 
stable quark matter with finite strangeness is only possible if the in-medium strange quark 
mass m* is much lower than the value we obtained above. The easiest way to achieve 
this is to choose a lower value of the current mass Wq. If we leave all other parameters 
unchanged, ttiq should be at most 85 MeV if we require n s ^ at the minimum of E/A. 
With this value we obtain much too low masses for K and rj (rriK — 390 MeV, m v ~ 
420 MeV), while E/A is still relatively large (1075 MeV). If we want to come down to 
E/A ~ 930 MeV we have to choose m s Q = 10 MeV or, alternatively, m s Q = 25 MeV and 
mJJ = Wq = 0. This is of course completely out of range. 
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We could also try to lower m* by choosing a smaller coupling constant G. (Since 
(uu){dd) is already very small at the densities of interest, m* is almost insensitive to K). 
However, with a lower G the vacuum masses of all quarks drop and correspondingly the 
vacuum bag pressure Bq. This causes the minimum of E/A to move to lower densities. 
It turns out that the corresponding chemical potential /z* drops even faster than m* and 
it is thus not possible to obtain stable SQM in that way. To avoid the strong density 
decrease we could increase the coupling constant K while decreasing G, e.g. in such a 
way that the vacuum masses of u- and d-quarks are kept constant. In order to get n s ^ 
at the minimum of E/A we have to lower GA 2 to 1.5 and to increase KA 5 to 21.27, almost 
twice the value of our original parameter set. For these parameters the energy per baryon 
number is still 1077 MeV. On the other hand we cannot further decrease the ratio G/K 
because that would flip the sign of the effective qq coupling in the pseudoscalar-flavor 
singlet channel, which is dominated by the combination 2G + ^K((uu) + (dd) + (ss)). In 
that case there would be no solution for the //-meson in vacuum. Hence there seems to 
be no realistic way to find absolutely stable SQM within our model. 

We could ask whether this result can change if we introduce additional interaction 
terms to the Lagrangian eq. (|l|). In particular vector interactions are known to be im- 
portant in dense matter. However, since vector mean fields are repulsive, the energy per 
baryon number will be even larger than before and SQM remains strongly disfavored com- 



pared with ordinary nuclei. The same is true for the medium effects described in ref. [p8 



which also tend to increase the energy. In fact, these effects are complementary to ours 
because they are most important (and most reliable) at high densities, whereas our dy- 
namical quark masses are most important at low densities. In ref. |T9 the chiral phase 



transition was studied within a "scaled" NJL model. There is no six-point interaction in 
this model, but the different flavors are coupled through a dilaton field. The coupling was 
found to be weak except for rather low choices of the gluon condensate P0| . This makes 
it very unlikely to find absolutely stable SQM within that model, although a quantitative 
study of this question has not yet been done. 

In summary, we investigated strange and non-strange quark matter within the NJL 
model. In the chiral limit the model behaves very similarly to the MIT bag model which 
has been mostly used in the literature to describe SQM. For realistic model parameters 
including finite current quark masses the NJL model shows a more complex behavior. This 
is due to the fact that bag pressure and effective quark masses are dynamical properties 
of the model and are therefore in general density dependent. We find that the dynamical 
mass of the strange quark stays very large at the relevant densities. As a consequence 
SQM is not favored compared with NSQM in the model and there seems to be no chance 
to find SQM with energies per baryon number lower than in ordinary nuclei. This almost 
rules out the original idea of absolutely stable SQM as a simple consequence of the Pauli 
exclusion principle. Of course we cannot completely exclude that more sophisticated 



mechanisms, like color superconductivity [29, 30], in particular the recently discovered 



mechanism of color-flavor locking pi| , change our results. This is, however, beyond the 
scope of this letter. 
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Figure 1: Properties of isospin-symmetric (n u = rid) stable quark matter as a function of 
r s = n s /riB in the chiral limit (parameters of ref. J2l|, but with m l = 0). (a) Energy per 
baryon number. The result of the present model (solid line) is compared with a bag model 
calculation with bag constant Bbm = 57.3 MeV/ fm 3 (dotted), (b) Dynamical masses of 
the strange (solid) and non-strange quarks (dotted), (c) Effective bag constant. 
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Figure 2: The same as fig. 1 but with finite current masses (parameters of ref. 
In panel (a) the result of the present model (solid line) is compared with bag model 
calculations with bag constant B BM = 105.2 MeV/fm 3 and m u = m d = 5.5 MeV, m s = 
140.7 MeV (dotted) or m u = m d = 52.6 MeV, m s = 464.4 MeV (dashed-dotted). 
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Figure 3: Properties of charge-neutral quark matter with electrons as a function of baryon 
number density ps (parameters of ref. f2lf). (a) Energy per baryon number in chemical 
equilibrium = p d = p u + /x e , solid line) and for non-strange matter (/i s = 0,p d = 
fi u + /i e , dotted), (b) Dynamical quark masses m u (dashed-dotted), m d (dashed) and 
m s (solid). The dotted line denotes the chemical potential /i s = \i d . (c) r\ = rii/riB for 
the different particle species: up quarks (dashed-dotted), down quarks (dashed), strange 
quarks (solid), and electrons (dotted). 



